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ABSTRACT 


This study explores the dynamics of diffusive 
convection which occurs when cold, fresh water overlies warm 
and salty water. Diffusive convection is generally observed 
in high-latitude regions, particularly in the Arctic Ocean, 
between the base of the mixed layer and the top of the 
Atlantic water. The primary convective regime in this 
region is characterized by the spontaneous formation of well 
mixed layers separated by thin high-gradient interfaces 
known as thermohaline staircases. 

Data analysis and analytical considerations are used to 
estimate the vertical heat/salt mixing rates and their 
dependencies on the large-scale environmental parameters. 
The ice-tethered profiler data from the Beaufort Gyre are 
analyzed to determine the origin of the thermohaline 
staircases and the mechanism for selection of the preferred 
layer thickness. Based on the analysis of Beaufort 
staircases, we suggest that the layer thickness, as well as 
the vertical heat/salt fluxes, is controlled by the patterns 
of merging events in which relatively small steps are 
systematically eliminated. Significant concerns are raised 
with regard to the direct extrapolation of laboratory 
derived flux laws to ocean conditions — a conventional 
approach to the analysis of the Arctic staircases. An 
alternative and more accurate method of analysis is proposed 
which involves recalibration of the laboratory-derived flux 
laws for the oceanic conditions. 

Extrapolated diffusive convective fluxes are in the 
range of 1-6 Wm''-2, which is comparable to magnitude of 


V 



fluxes currently unaccounted for in the Arctic heat budget. 
We propose that the parameterizations of the diffusive 
fluxes in thermohaline staircases, developed in this thesis, 
can be used to enhance understanding of Arctic climate 
changes and predictive capabilities of large-scale numerical 
models. Preliminary findings are indicative of the 
importance of diffusive convection for sound propagation in 
the Arctic region — the problem of great interest for 
various Naval research applications in the area. 
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I. 


INTRODUCTION 


The polar ice cap, consistently regarded as a hindrance 
to data collection in the Arctic, is currently under¬ 
sampled. With evidence of increasing Arctic temperatures 
and decreasing sea-ice thickness, the need to parameterize 
unresolved ocean processes, and include them in climate 
models, is becoming a necessity. This study focuses on the 
double diffusive region of the Beaufort Gyre. 

A. DOUBLE DIFFUSION 

Before discussing the Arctic data, the basic concepts 
behind double diffusion should be reviewed. Although double 
diffusion occurs in many different systems, it has a vast 
impact on ocean dynamics. Turner (1973) provided an entire 
chapter on the dynamics of double diffusive regimes in his 
book on the buoyancy effects in fluids. He describes 
double-diffusion as the process related to the presence of 
more than one component contributing to the density of the 
ocean with different molecular diffusivities. In the ocean 
case, the dominant density components are temperature and 
salinity. Ocean double diffusion takes two forms. The most 
studied is the salt finger regime, where warm, salty water 
lies above cold, fresh water. The other regime, which is 
the focus of this study, is diffusive convection where cold, 
fresh water overlies warm, salty water. Stratification 
patterns which support salt-fingering and diffusive 
convection are shown in the schematic diagrams in Figure 1. 
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Figure 1. Simplified diagram of ocean double diffusive 

regimes, a) Salt Finger Regime b) Diffusive Convection 

Regime. 


1. Salt Fingers 

The majority of research has focused on the salt¬ 
fingering regime. Stommel et al. (1956) brought interest to 
the salt finger regime by hypothesizing that by exploiting 
the warm, salty water over cooler, fresher water, an energy 
source termed "the perpetual salt fountain" could be formed. 
The concept being that a thin membrane in the shape of a 
tube could allow temperature diffusion, but not salinity 
diffusion, and create an energy source by making a parcel of 
salty water rise as it tries to reach a state of 
equilibrium. While the proposed idea has yet to be fully 
developed in the ocean, the interest in salt fingering made 
its debut. 

Stern (1960) described salt-fingers as instabilities 
that arise when relatively warm, salty water lies over 
colder, fresher water of a higher density and consists of 
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long fingers of rising and sinking water. The warmer, salty 
water rapidly loses heat while retaining its salinity. The 
increasing parcel density makes it sink further, drawing 
more warm, salty water from above, resulting in the fingers 
of salty fluid. Figure 2 shows the resulting salt fingers. 
In this case the magenta dye is applied to the warm, salty 
water. 



Figure 2. Laboratory simulation of salt fingers. The 

magenta plumes are salty water pushing into the colder 
fresher water. (From University of Alaska Fairbanks 
website, www.ims.uaf.edu). 

2. Diffusive Convection 

The dynamics of double-diffusive convection differ 
substantially from those of salt fingering. Diffusive 
convection occurs in both a smooth gradient region and a 
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stepped configuration. 


which is more common in nature 


(Figure 3) . For the diffusive convection case, relatively 
cold, fresh water overlies warm, salty water with a thin 
interface in between. 


Cold, Fresh 

Thin interface 


ir Ah 

Warm, Salt 

/ 


Figure 3. Basic set-up of diffusive layers. 


Molecular fluxes across the interface are: 




aAT 

Ah 


and = k. 


PAS 

Ah 


where a is the coefficient 


of volume expansion due to unit temperature change, P is 
the proportional density change produced by unit density 
salinity change, Fj.j^ and Fg„ are the molecular fluxes of heat 

and salinity respectively, and /fj. and are the molecular 
diffusivities of heat and salt, respectively. 


Since , for aAT and PAS of the same order of 

magnitude, the density flux Fp= pPP^^ -paFj^^ <Q. The rapid 

molecule diffusion of temperature, relative to salt, 
produces a downward density flux across the interface. The 
region immediately below (above) the interface becomes 
denser (lighter), which maintains the top-heavy convection 
in mixed layers. Figure 4 illustrates the physical 
mechanism for diffusive layering. 
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Figure 4 . 


Physical mechanism behind diffusive layering. 


There are several regions of the world's oceans with 
the potential for double-diffusive convection, including the 
Arctic Ocean (Kelley et al. , 2003) . The intensity of 
diffusive convection is controlled by the so-called density 
ratio, which, for the diffusive convection case, is defined 
as: 


p^s 

«Ar 


( 1 ) 


where: 


a = - 

P 


dp 

ydT* 




P 


dS 


\ 

J S*,P* 




( 2 ) 


Fluxes decrease with increasing density ratio. The 
range is most susceptible to diffusive convection 

(Kelley et al. , 2003). Analysis of Arctic density ratios 
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shows that the majority of the Arctic Ocean, specifically 
the Beaufort Gyre, has density ratios in the 3-10 range. 
(Figure 5) 



Figure 5. Regions susceptible to double diffusion. The 

light grey areas indicate where density ratio is from 
3-10. Dark grey indicates l<i?^<3. The red circle 
indicates the Beaufort Gyre. 

The Arctic is diffusively stratified due to the process 
of ice formation, brine rejection, and subsequent melting. 
The process creates the cold, fresh layer that lies on top 
of the warm, salty water carried from the Atlantic Ocean to 
the Arctic, through typical ocean circulation. Figure 6 
illustrates Arctic currents and circulation patterns. 
Figure 7 provides a vertical cross section of the Arctic 
Ocean. 
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Figure 

(I 


6. Arctic circulation pattern, 

conditions favorable for diffusive 
llustration by Jack Cook, Woods Hole 

Institution.) 


resulting in 
convection. 
Oceanographic 
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Wlwly fofined ice. 

Saltis I 
released 


\ 


SO meters 


Denser, 
salty water 
sinks 


Cold, 

fresh 

water 


Cold, 

saltier 

water 


200 meters 


Warm, salty 
Atlantic water 


Figure 7. Cross section of Arctic water mass. Cold, 

fresh water overlies warm, salty Atlantic water. 

(Illustration by Jayne Doucette, Woods Hole 
Oceanographic Institute.) 

In Figure 8 Fernando (1987) simulated the layering 
process by heating a salinity gradient from below. The 
diffusive layers are indicated in the diagram. Similar 
laboratory simulations have also been conducted using a 
sugar solution as the denser liquid. Figure 9 is a more 
completely evolved set of diffusive layers in the 
laboratory. Turner (1968) used this shadowgraph to show the 
scale of layers. 
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Figure 8. Side view of diffusive layering created by 

heating a salt gradient from below, after Figure 6 
(Kelley 2003) A is the first quasi-stationary 
interface. B is the secondary interface, forming above 

it. 



Figure 9. A series of convecting layers and heat-salt 

diffusive interfaces (Turner 1968). 
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B. THEORETICAL BACKGROUND 

Turner (1965) hypothesized that if steps are 
sufficiently large {H^cc), the dependence of the vertical 
fluxes of temperature and salinity (iy and iy) on H would 
be weak, and therefore fluxes would be determined by the 
temperature and salinity variations (AT,AS) across the step: 

Fj. = (aAT, PAS,k^ ,K^,v,g) = F^ (aAT, R^,/Cj.,g,T,PJ 
F^ = Fg (aAT, PAS,,/c^,u, g) = Fj. (aAT, R^,/Cj.,g,T,Pp ^ ^ ^ 


where u is viscosity, r is the diffusivity ratio and is 
the Prandtl number 




J 


Using the dimensional argument that a non-dimensional 
number can depend only on other non-dimensional numbers. 
Turner then suggested that the so-called Nusselt number. 


Nu = 


Ft 

AT 


measuring the fluxes would be related to the 


Raleigh number, measuring boundary forcing: 


gaATH^ 

KjO 


Nu = C(RJ^aF^ 


^ KjV ) H 


(4) 


From (3), we know that there should be no dependence on 

4 

H in (4), therefore n=l/3 or aFj~(aAT)^. These equations 
are called the "4/3 Flux Laws". (Turner 1973) 
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Another non-dimensional quantity of interest is the 




flux ratio: y = —which measures the fraction of the 


potential energy lost by temperature stratification that 
goes into the potential energy of salt stratification. 


Turner (1965), argued, on the basis of dimensional 
considerations and laboratory experimentation, that for 
given molecular properties, the flux ratio is a function of 
Rp . Experiments yield a striking result that for Rp>2 the 
flux ratio remains constant to within the experimental 

BF 

error. In the laboratory this mean value is -^ = 0 . 15 , with 

aFj. 

maximum heat density flux for a given AT occurring at = 2 . 

Figure 10 is a plot buoyancy flux ratio versus density ratio 
in the laboratory. 
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Figure 10. Buoyancy flux ratio Rp(R^) for diffusive 

convection in the laboratory. Data from three 
different laboratory sources Crapper (1975), Turner 
(1965), and Newell (1984) from Kelley (1990). The 
solid line is the empirical fit. 

The aforementioned models offer simple conceptual 
explanations for the dependences of vertical fluxes on 
temperature and salinity variations across the steps 
(AT,AS) . It follows that knowledge of what determines AT 
in the staircases is necessary. Kelley (1988) and Radko 
(2005) suggest that a layer-merging process controls step 
characteristics. This merging process, which will be 
discussed further in Chapter III, can be defined as the 
coalescence of thin, adjacent layers to form new, thicker 
layers, and their eventual equilibration. The major goal of 
this study is to examine the validity of flux law models and 
predictions of merging theory. For that, we turn to Arctic 
data. 
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c. 


ICE-TETHERED PROFILER 


Originally, thermohaline staircases were not visible in 
ocean profiles. This was due to data taken only at standard 
levels versus a continuous profile of the region. In order 
to characterize the diffusive convective layers, high- 
resolution data are necessary. This study uses high 
resolution data collected from the Woods Hole Ice-Tethered 
Profiler Project. 

In conjunction with the Beaufort Gyre Exploration 
Project, a team at Woods Hole Oceanographic Institute (WHOI) 
has launched several Ice-Tethered Profilers (ITP's) in the 
Beaufort Sea. Analysis of diffusive convection primarily 
uses data from ITPl and ITP 3. 

ITP's are moored profilers that are designed to make 
daily measurements of water, pressure, salinity 
(conductivity), and temperature in a region where ice 
impedes access and collection of this data. Designed to 
last for three years when deployed, as long as the ice floe 
survives, the profiler sends its measurements and GPS 
location to a surface float, which is then relayed back to 
WHOI. The profiler is made up of a small surface capsule, an 
800 meter line, a moored profiler, and a small anchor. The 
profiler is modified so to fit through a 10" hole in the 
surface (Krishfield et al., 2006). 
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Figure 11. Schematic of Ice-Tethered Profiler 

(Krishfield et al., 2006). 

The high resolution (25 centimeter) data posts to the 
WHOI ITP website for further analysis in its raw form. 
There are future plans to de-spike, filter, and correct the 
data using a sensor response correction (Johnson, 2007) . 
These corrections should occur in the fall of 2007 and be 
considered for future work. 

ITP 1 was the second of seven ITP's deployed. It was 
launched in the Beaufort Sea on August 15, 2005 at 78° 51.1 
N, 150° 15.9 W during the JWACS 2005 cruise on the CCGS 
Louis S. St. Laurent. ITPl was deployed on a 4.6 meter thick 
multiyear ice floe in a location designed to drift through 
the Beaufort Gyre Observing System (BGOS) array over the 
following year. ITP3, on August 23, 2005 at 77° 36.1 N, 
142° 11.8 W, during the same cruise. It was deployed on a 
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3.5 meter thick multiyear ice floe in a location designed to 
drift through the Beaufort Gyre Freshwater Experiment (BGFE) 
array during the following year. Both units operate on a 
fast sample schedule of 4 one-way profiles between 10 and 
760 meters depth each day. 



[ 5000 
4000 
3000 
2000 

r 

1000 
500 
200 
•-L 0 



ITP drift (yellow line) and latest location (triangle), 
BGOS moorings (white circles) and annual Ice drift from 
lABP (grey vectors) on IBCAO bathymetry (shading). 


Figure 12. ITP 1 and ITP 3 drift tracks, (a) Drift track 

for ITP 1 and (b) is ITP3 drift track. From Woods Hole 
Oceanographic Institute Ice Tethered Profiler site 
(http://WWW.who!.edu/itp/data.html). 
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II. CHARACTERISTICS OF THERMOHALINE STAIRCASES 


A. STEADY CHARACTERISTICS 

The temperature, pressure, and conductivity raw data 
were used to create temperature-depth and salinity-depth 
profiles through the diffusive layer, characterized by 
prominent staircases. In the Beaufort Gyre Region, 
"staircase" structures were prominent from depths of about 
200 to 400 meters. Therefore, all analyses were computed 
within that depth range. Figures 13 and 14 are temperature- 
depth and salinity-depth profiles illustrating the staircase 
characteristics. As denoted in Figure 13, H is defined as 
the distance between the midpoints of adjacent interfaces. 
Figure 14 indicates A5 measurements. AT is found in the 
same manner. All vertical calculations were conducted with 
this data. The step is defined as the vertical rise between 
two adjacent interfaces. a and [5 were calculated using 
the MATLAB seawater toolbox. Although they vary with depth, 
average values of a are O^IO^'-S and P are 8*10''-4 

(PSU)”^. Average AT was 0.04 degrees Celsius, and average 

A5 values were 0.01 PSU. The density ratios R = 

[ aAT 

computed for the ITP data, ranged from 3-7. Table 1 shows 
the data for ITP 1 and ITP 3 averaged by month. 
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Figure 13. Temperature Depth Profile. The section of 

the profile is highlighted. H calculations are 
determined from the midpoint of each staircase depth. 
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Salinity (psu) 


Figure 14. Salinity depth profile. Same location as 

previous figure. (Salinity is calculated using MATLAB 
polyfit to de-spike data). 
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AT 

H(m) 

Rp 

ITP1 




September 

2005 

0.0403 

4.1635 

3.2713 

October 


0.0421 

4.7896 

3.6998 

November 


0.0437 

4.6067 

3.8795 

December 


0.0466 

4.2825 

4.759 

January 

2006 

0.0437 

4.1593 

3.8973 

February 


0.0434 

3.8276 

3.8536 

March 


0.0456 

4.3723 

3.7838 

April 


0.0406 

4.0872 

3.8512 

May 


0.0397 

3.4183 

4.0875 

June 


0.0429 

3.9029 

4.1846 

July 


0.0472 

4.31 

4.3198 

August 


0.0441 

3.9876 

4.6595 

September 


0.0435 

4.0714 

4.9589 

October 


0.0466 

4.2825 

4.759 

ITP3 




September 

2005 

0.0373 

3.2653 

5.5181 

October 


0.0432 

3.6059 

4.548 

November 


0.0394 

3.5775 

4.6802 

December 


0.0394 

3.4346 

4.8329 

January 

2006 

0.0378 

3.0995 

4.7914 

February 


0.04115 

3.3285 

4.638 

March 


0.0455 

3.8186 

4.4491 

April 


0.046 

3.5688 

4.2637 

May 


0.0455 

3.6922 

4.3631 

June 


0.0477 

4.1085 

4.734 

July 


0.0356 

2.7624 

5.4397 

August 


0.0391 

3.1684 

5.3817 

Average over both profilers 

0.0426 

3.8342 

4.446 


Table 1. 


Average data 
over 


from ITP 1 
the entire 


and ITP 3 by month and 
data set. 
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Figure 15 plots aAT (non-dimensional temperature 
variation with depth) as a function of H (calculated step 
height.) Although only the month of January is presented, 
similar results were obtained for each month. From this 
comparison, it is theorized that the change in temperature 
is relatively linear, and is nearly independent of the 
height of the step. The slightly negative slope is caused 
by the dependence by a on pressure. In general, H values 
increase from about 2 meters at a depth of 240 meters to 
above 20 meters at a depth of 400 meters, while AT remains 
on average of 0.04 degrees C. This observation is 
consistent with Turner's suggestion that the vertical fluxes 
of heat and salt are independent of step height (Equation 
(3) ) . Indeed, if the vertical fluxes Fj and are 

vertically homogeneous, then AT, AS = constant. Figure 16 
shows a similar relationship between salinity and H. 


X 10 * January 2006 (Potential Temp ITP1) 



Figure 15. Non dimensional aAT is plotted as a function 

of layer thickness, H (meters). 
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Figure 16. 


Non-dimensional yS’AS' is plotted as 
of layer thickness H (meters). 




function 
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aAT 


In Figures 17 and 18, alsT and ySAS are presented as a 


function of density ratio, . 


X 10® 



Figure 17. Comparison of aAT and using average data 

from September 2005 to October 2006. 
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3.2 3.4 3.G 3.B 4 4.2 4.4 4.G 4.8 5 


Density Ratio 

Figure 18. Comparison of /51S.S and using average data 

from September 2005 to October 2006. 

In Figure 19, the density ratio over the interface is 
compared with the density ratio over the step. For this 
comparison, density ratio is first calculated over the step, 
and then calculated using two points on the adjacent 
interface. The density ratios for the interface and steps 
are then averaged by month and plotted. The increased value 
of ^pinterface oveu Rpstep^^ ou Indlcatlon of the ambiguity in the 

definition of density ratio, and leads to a possible source 
of error in the 4/3 flux laws. 
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ITP1 data averaged by month from September 2005 to October 2006 



Figure 19. Comparison of density ratios over the 

interface and density ratios over the step. 


There is general agreement in the literature that the 
double-diffusive temperature and salinity fluxes are 
controlled by the density ratio. While this relationship is 
known to exist, there is still considerable uncertainty in 
both the dynamics and values of fluxes in the thermohaline 
staircases using a straightforward extrapolation of the lab- 
based flux estimation. Padman and Dillon (1987) calculated 
thermohaline fluxes in the Beaufort Gyre to be about 0.1 
Wm'^-2 for 2<R^<1. A comparison between the thermohaline 

staircases that Padman and Dillon explored and the current 
data indicate that in the 20 years between data collection, 
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there were significant changes in the water mass. This 
leads to changes in the calculated fluxes. 

In 1987, the depth range of the staircases was 300-430 
meters with AT on order of 0.01 degrees C vice the 0.04 
found today. There is a more fundamental difficulty with 
regards to the earlier estimates of fluxes in diffusive 
staircases. In addition to the differences in water column 
characteristics, Padman and Dillon's calculations were made 
using direct extrapolation from laboratory numbers. In the 
salt finger regime, such an extrapolation was long ago 
proven unreliable. In the diffusive convection regime, 
Kelley (1990) also suggested that direct application of the 
4/3 Flux Law to ocean conditions is equally unreliable. In 
order to validate this theory, Radko (2007)'s merging 
theorem was applied to two components, temperature and 
salinity. These theoretical values, using laboratory 
results, were compared to the actual data. The result was 
confirmation that the 4/3 Flux laws cannot be directly 
extrapolated from the laboratory setup to diffusive 
convection in the ocean. 
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B. 


EVOLUTIONARY PATTERN 


Initially, the possibility of mergers existing in the 
ITP data came from a plot of conductivity versus time. 
Conductivity was used vice salinity because it could be 
taken in its raw form. In Figure 20, lines of constant 
conductivity, approximating the observations, are 
superimposed on the measured data. The merging events 
occurred twice in this figure — that is where the lines of 
constant conductivity meet and create fewer lines. 


ripi 



200 4X 600 800 1000 1200 1400 

Time 


Figure 20. Examples of merging events, as seen from the 

inspection of the conductivity data. 
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Overall, we identified ten clear mergers to use in the 
following (Chapter III) theory-based analysis. Several of 
these mergers are presented in Figure 21. Visual inspection 
of all merging events suggest that layers merge when 
slightly stronger, in terms of temperature and salinity 
variations, interfaces grow, while weak interfaces decay and 
ultimately disappear. This is the evolutionary pattern 
classified by Radko (2007) as the B-merger scenario. 


** , . 





1 . .. ............ : 
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Figure 21. Temperature variation in time. Each data 

point represents temperature on a particular profile. 
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III. MERGING THEOREM FOR A TWO COMPONENT FLUID 


Radko (2007) theorized that the formation of diffusive 
layers starts with multiple, small-scale, homogenous density 
layers that merge until they reach equilibration. There are 
two types of merging processes. Radko defines them as H- 
merger and B-merger. The H-merger is defined when slightly 
thicker layers grow while thin layers shrink and ultimately 
disappear (Figure 22). In the case of B-merger (Figure 23), 
stronger interfaces grow at the expense of weaker 
interfaces. While theoretically both types are possible, 
data indicates that only the B-merger exists in the ocean, 
therefore the included derivations are limited to the B- 
merger. 


H-Merger 




t=3 


Figure 22 


Schematic of H-merger evolution in time, 
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B-Merger 



1=2 


t=3 


Figure 23. 


Schematic of B-merger evolution in time 


The derivation for one component and one-dimensional 
mergers in Radko (2007) was modified to two components, 
salinity and temperature, yielding the following description 
of the B-merger events. This is accomplished by conducting 
a linear stability analysis for a series of identical steps 
and using it to formulate general criteria for B-merging 
events. 

In order to establish the basic principals of the 
staircases evolution in time, we consider the stepped 
solution of the one-dimensional conservation equations: 

dT _ dF ^' 

dt dZ , n > 

. r 

dS _ dF^ 

dt dz . 


The stepped solution is shown in Figure 24. Figure 24a 
represents a staircase consisting of identical steps, the 
basic steady-state. Figure 24b represents a slightly 
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perturbed state in which steps are no longer equal. The 
temperature-salinity variations across the two steps in 
Figure 24(b) are expressed as: 


T 

^\2 

T 

-*23 


-T,,=T-d 
-Tn=T + S 
-S,,=S-s 
— iSj2 — S s, 


(S,s)« (T,S) 


where T = (dT /dz)H and S = (dS /dz)H are the variations over 
step in the undisturbed staircase in Figure 24 (a) and (t 
are small perturbations. Assume that H is constant 
equal where H is the distance between the interfaces 
z = z„ and z = z„^i . 


( 6 ) 


one 

.S) 

and 

at 
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Figure 24. Schematic diagram illustrating the stability 

analysis for an infinite series of interfaces. (a) 
Basic state consisting of identical steps, (b) 
Perturbed state in which the T,S jumps at the even 
interfaces are slightly decreased, and the jumps at odd 
interfaces are increased. 


Then integrating (5) over the interval [zbot^ Ztop] where 
Ztop and Zbot are two levels which do not vary in time: 



^bot ^bot 
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Using the identity: 



( 8 ) 


Equation (7) is rewritten as 

^hoi 

> 

rl 

-I Sdz=F‘(z„,)-FUz^,) 


Applying the integral relation in equation (9) to 
individual layers [z„,z„+i] yields, where Fi=F(zi): 



( 10 ) 


The major contribution to the integral of buoyancy in 
(10) comes from the interior part of layers where the 
buoyancy gradient is nearly uniform, and therefore. 


71+1 

-«+l 

I Sdz^S^„,H 


^«+l 


> (11) where and are the 


temperature and salinity values at the layer centers. 
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Applying (8) and (9) to two consecutive layers 
where n=l, 2 results in: 

Using equations (6) and (12), the equations for 
evolution of d and s can be formulated, and then simplified 
using periodicity conditions. 



at 

H^ = -1{F^-F,^) 
at 


(13) 


Considering weak perturbations such that (S ,s H) and 
linearizing (15) yields: 



(14) 
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At this point, it becomes necessary to specify how the 
fluxes in (14) depend on temperature-salinity jumps across 
the steps. For that, we replace the instantaneous fluxes in 
interfaces iy in (15) with the steady 'one-step' flux 

F (Radko 2007) based on the parameters of a region extending 
between the centers of two adjacent layers andz^^ j) : 


p! 

Pi 

Pi 

Pi 


~P (^23 “^ 12’*^23 

«F^(r,2-r„„5,2 

«#^(r,2-r„,,5,2 


5'oi, Zi2 
*^12’ -^23 
*^ 01 ’ -^12 
*^01 ’ -^12 


-^oi) 

-^ 12 ) 

-^Ol) 
~ ■^01 ) 


(15) 


The difference between fluxes at the adjacent 
interfaces in (14) reduces to: 


F^ 

F^ 
^ 2 


Fj =15 


F/ = 15 


dF^ 

df 

dP^ 

df 


+ ls 


+ ls 


dF^ 

~w 

dP^ 

dS 


and (14) becomes: 

d 45 dP^ 4g dP^' 

dt ^ H df H dS 

d __45 dP\ 4s dP^ 

dt^^ H df H dS . 


(16) 


(17) 


Substitution of the normal modes ((5’,^,/z) = ((5g,£'g,/zQ)exp(/l0 

in (17) yields the eigenvalue equation for the B-growth 
rates: 


t+- 

"" H 


dP^ dF 

+ 




dT dS 


2 16 

B ^2 


dF^ dF^ 
df dS 


dP^ dP^ ^ 

df dS , 


(18) 
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If the free coefficient is small (to be expected since 
it is proportional to the flux ratio variation) the positive 
root can be further approximated by: 

dP^ dP^ dP^ dP^ 

"" H 

df dS 

The tildes denote the equilibrium one-step parameters. 
The latter affords an attractive opportunity to validate the 
merging theory (Huppert, 1971; Kelley, 1988; Radko, 2005) by 
comparing the growth rates of the observed merging events 
with the theoretical prediction for F{R^) . 
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IV. APPLICATIONS OF THE MERGING THEOREM TO OCEANIC 

DATA 


In order to apply the merging theorem (19) to actual 
field data, explicit expressions for fluxes are needed. 
While little is known about the flux laws in the ocean, we 
assume that the general structure of the laboratory-derived 
flux laws pertain to the ocean, but additional calibration 
of amplitude is necessary. Adopting the derived flux laws, 
but recalibrating them for oceanic conditions, (e.g. Kelley 
1990) yields: 


aF' =AC 


( 7 2V'^ 

V ^ y 


/3F^ =aF^y, / = 


iaATy'^\ C = 0.0032 exp(4.8/i?°"') 


( 20 ) 


"A" is an adjustable coefficient included to 
recalibrate Kelley-Turner flux laws for the ocean. 


An alternative lab-based formulation was suggested by 
Marmorino and Caldwell (1976): 


= {Qm5)kp 


ga 

yKv) 


AT 


4/3 


( 21 ) 


H 


H 


— = 0.101exp|4.6exp[-0.54(i?^-l)' 

SP 


(22: 


Finally, substitute (20) and (22) into (19) and obtain 
specific solutions. Simplifying (19) takes the form: 




( 23 : 
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Now, we apply the ITP data. The exponential variation 
of AT,AS across the merging layers confirms the 
interpretation of the merging as an instability process. 
The growth rate of each merger is determined from the best 
fit of the exponent exp(/lO to (AT^-AT^i^), where AT^, is the 
temperature variation in the merging layer, and AT^^ is its 
value at the beginning of the merging event. Two examples of 
growth rate determination are provided. In the first 
example, the density ratio is slightly higher, leading to a 
smaller growth rate and a longer period for merging (Figure 
25) . In the second example, the density ratio is only 
slightly smaller, but the growth rate is twice as large, 
yielding a much shorter merging period (Figure 26) . 



Time in Days 

Figure 25. Change in temperature vs. time during the 

merging period. The normalized growth rate for this 
merger is 4.7E-04 and the density ratio is 2.8631. 


38 












Figure 26. Change in temperature vs. during the merging 

period. The normalized growth rate for this merger is 
8.6E-04 and the density ratio is 2.4035. 


The merging growth rates are then normalized as 
follows: 


1 . 


XH 


'normal 


( 24 : 


where H is the initial step height and oc = 6e — 5K _ They are 
then compared with the normalized growth rates from the 
theoretical/laboratory values. Kelley (1990) suggested that 
the coefficient C in the 4/3 flux law differed based on 
density ratio. Based on laboratory estimates. 


C =.0032exp 


^ 4.8 ^ 


pO.72 


{25] 


This differs 
suggested 


from Marmorino and Caldwell 

C=0.085*0.101*exp(4.6exp[(-0.54{R^-l}]) . 
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The differences in the calculations of laboratory data 
lead to questions regarding the precision of laboratory 
measurements. In order to estimate the error of the model 
predictions, both Kelley and Marmorino and Caldwell's 
formulations were independently used. Figure 27 is based on 
Kelley's expression, which required a coefficient A=30.008 
to fit Beaufort Gyre data. Marmorino and Caldwell's 
coefficient, was A=15.75 and the corresponding prediction 
for is shown in Figure 28. Table 2 shows a numerical 
comparison of the fit. 


Data vs. Theory 



Figure 27. Using the coefficient A=30, theoretical data 

fits Beaufort Gyre data within 95% prediction bounds. 
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X 10^ merging theory 



Figure 28. Using the coefficient A=15.75 applied to 

Marmorino and Caldwell (1976), theoretical data fits 
Beaufort Gyre data within 95% prediction bounds. 


Rp 

AT 

H 

A normalized 

A normalized 
A=30 

A normalized 
A=15.75 

1.9106 

0.1572 

24 

2.70E-03 

2.58E-03 

2.48E-03 

2.2509 

0.1007 

7 

3.1 IE-04 

3.1 IE-04 

7.93E-04 

2.2737 

0.087 

3 

3.60E-04 

3.60E-04 

7.43E-04 

2.4031 

9.61 E-02 

4.6 

1.06E-03 

7.17E-04 

5.24E-04 

2.4035 

1.12E-01 

3 

8.63E-04 

7.17E-04 

5.23E-04 

2.4267 

1.18E-01 

5 

6.08E-04 

6.83E-04 

4.93E-04 

2.8631 

6.99 E-02 

4 

4.71 E-04 

3.16E-04 

1.90E-04 

3.2997 

8.20E-02 

3.7 

1.28E-04 

1.75E-04 

9.25E-05 

3.3728 

7.91 E-02 

2.8 

2.81 E-04 

1.60E-04 

8.35E-05 

3.3879 

9.09 E-02 

3 

1.90E-04 

1.58E-04 

8.18E-05 


Table 2. Data comparison of calculated normalized growth 
rate from actual data and the normalized growth rates 
using a constant factor with the 4/3 flux law. 
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The coefficient A, which we believe is necessary to 
convert from laboratory simulations to ocean data leads to 
relatively high heat fluxes low density ratios 

(1.9 < < 3.4 ) . When applying the A-coef f icient to 

characteristics of multiple thermohaline staircases, the 
average flux is 6.65 Wm''-2 for Kelley's formulation and 
3.717Wm''-2 for Marmorino and Caldwell's. The range of 
fluxes is obtained using ITP 1 data from April of 2006. 
Average AT, , and H values were calculated in 20 meter 

increments from 220 to 360 meters (except for 320-360 meters 
where steps were too large to make accurate calculations 
from a 20 meter interval). This range covers 1211 steps. 
Table 3 shows the data averaged by layer. 


Depth in meters 

220-240 

240-260 

260-280 

280-300 

300-320 

320-360 

AT 

0.0426 

0.0405 

0.0423 

0.0635 

0.0439 

0.0368 

AS 

0.0132 

0.0125 

0.0118 

0.0202 

0.0145 

0.0138 

Rp 

6.1626 

5.5163 

4.6419 

4.9691 

4.9435 

5.6985 

H (meters) 

2.3520 

2.6326 

3.1538 

4.4798 

6.0994 

5.9379 

Total Steps 

317 

306 

275 

157 

70 

86 








Kelley's Flux Law 







HFw/A=1 

0.1570 

0.1637 

0.2134 

0.3327 

0.2045 

0.0126 

HFw/A=19.6078 

3.0783 

3.2091 

4.1853 

6.5227 

4.0099 

0.2472 

HF w/A=30 

4.7090 

4.9099 

6.2648 

9.9797 

6.1351 

4.1831 

HFw/A=31.01522 

4.8692 

5.0761 

6.6201 

10.3175 

6.3427 

0.3910 








Marmorino and Caldwell’s 
Flux Law 







HFw/ A=1 

0.1524 

0.1606 

0.2176 

0.3370 

0.2076 

0.1366 

HFw/A=8.75 

1.3333 

1.4056 

1.9039 

2.9488 

1.8162 

1.1953 

HF A=15.75 

2.4000 

2.5300 

3.4270 

5.3079 

3.2691 

2.1515 

HFw/ A=17.955 

2.7360 

2.8842 

3.9068 

6.0510 

3.7268 

2.4527 


Table 3. ITP 1 data for April 2006 in 20-40 meter layer 

increments. 
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Table 4 uses the averaged data from Table 1 to provide 


a range of possible heat fluxes in the Arctic. The range 
varies significantly based on the size of the coefficient A. 
This suggests that further study is necessary to determine 
the validity of the 4/3 flux law when applying it to ocean 
conditions. 


Averages 
from Table 1 


Original 

HF Wm''-2 

Range of A 

HF Wm''-2 

Low Bound 
HF Wm''-2 

Best Fit 

HF Wm''-2 

High Bound 
HF Wm''-2 

AT 


Marmorino 

0.2360 

8.75<A<17.96 

2.0652 

3.7174 

4.2378 

H(m) 

3.8342 

Kelley 

0.2216 

19.61<A<45.90 

4.3459 

6.6492 

10.173255 

Rp 

4.446 








Table 4. Heat Fluxes in the Beaufort Gyre using the 

calibrated flux laws. 
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V. 


ACOUSTIC APPLICATIONS 


A. EFFECTS OF THERMOHALINE STEPS ON ACOUSTIC TRANSMISSION 

The "fine-structure" and "microstructure" 
characteristics of diffusive convection described above were 
rarely considered to have an impact on acoustic 
transmission. However, with recent high-resolution ocean 
data illustrating the depth and range that thermohaline 
staircases cover, these small-scaled phenomena are given new 
consideration. Chin-Bing et al. (1994) studied the effects 
of the thermohaline staircase step-structure on ocean sound 
velocity profiles (SVP's) and then performed acoustic 
propagation simulations using the oceanographic data. 

The simulations done by Chin-Bing et al. , used data 
collected from the northeast coast of South America in the 
fall (October-November) of 1985. The staircases studied 
were the result of salt-fingers and therefore had different 
characteristics (i.e., step size and density gradients) from 
the present studies of thermohaline staircases. However, 
these results led to the initial consideration of acoustic 
propagation through diffusive convective layers. 

B. INITIAL APPLICATION TO DIFFUSIVE CONVECTION 

Wide-angle capable, research-level. Parabolic Equation 
(PE) models were used to simulate long-range transmission 
loss through ocean thermohaline steps. The wide-angle models 
were used due to the deep ocean waveguide which supports 
fully retracted propagation for angles of ±15° . Due to the 
sensitivity of PE models to range step and depth step, the 
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model was run decreasing range and depth steps until the 
results converged on a stable answer. This led to the use 
of range and depth step sizes of 1/24 the wavelength. The 
frequency choice of 50Hz gave an acoustic wavelength on the 
same order as the step-size in the SVP. The source depth 
was 400 meters, placing it in the mid-depth of the steps. 
Results were plotted for receiver depths of 200, 400, and 
800 meters, that is, above, in, and below the step region 
respectively. In order to determine the importance of 
transmission loss through the staircases, a comparison was 
made between the transmission loss of the high-resolution 
step sound velocity profiles (SVP) and lower resolution data 
only taken at standard depths. 



Figure 29. a) Typical temperature and salinity plot for 

salt finger region. Salinity, Potential Temperature 
and Density plots show thermohaline step 
characteristics. b) Sound velocity profiles showing 
the step data, standard depth profiles and historical 
averaged sound velocity profile data. From Chin-Bing et 

al. (1994) . 
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The resulting simulation showed that when the receiver 
was above the step region there was very little difference 
when comparing the transmission loss of the two SVP data 
sets. The same determination was made when the receiver was 
under the staircase structure. There was, however, a 
significant difference within the step region--5 decibels 
over a 10-kilometer range. The greatest impact on the ray 
paths were when both the source and receiver were at 400 
meters, both were within the staircase structure. 

The numerical experiment then went on to test the 
effects of the staircases on backscatter. Using the same 
50Hz source, a source depth of 100 meters below the surface 
was chosen. Chin-Bing et al. found that in the low 
frequency regime of the study, the staircases tend to have a 
negligible effect; however, they felt as the frequency of 
the acoustic field increased so should the magnitude of the 
backscatter. They surmise that given frequencies in the 
sonar range i.e. 3.0 kHz, the backscatter from thermohaline 
steps could be significantly important. 

C. THE ARCTIC CASE 

From the above research, we thought it prudent to test 
the impact of the diffusive convection form of thermohaline 
staircases, given the probability of increased interest in 
the Arctic. The raw data from ITP profile 1013, chosen for 
its well formed staircase structure, was used for basic 
experimentation (Figure 30) . 


47 




Figure 30. Temperature and depth profile used for 

acoustic simulation. 

Raw temperature data and converted salinity (from 
conductivity) data was used for the idealized experiment. 
Using a basic fit of the data, the staircase structure was 
smoothed from the profile in order to offer a comparison in 
the effects of the microstructure on sound propagation. 
Using Del Grosso's Equation (Del Grosso 1974), sound speed 
profiles were calculated for simulation purposes. Figure 31 
shows the two sound velocity profiles. The only region of 
the profiles that was smoothed was the staircase region, so 
that effects of other factors were not included in the 
calculated impact. 
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SSP1013 / smaothed / with and without stanises 



Figure 31. Sound Velocity Profile from ITP 1 rawl013 

data. Red profile smoothes staircase structure, blue 
includes staircases. The graphic on the right is an 
enhancement to the SVP so that the staircase region is 

more visible. 

The parameters for the simulation were a planar, 
perfectly reflecting surface. (This is not the case for 
under the ice sheet, where there would be increased 
scattering of sound as the ray paths encountered the ice.) 
The bottom was simulated as a 17 meter thick sediment layer 
with a bottom half space. A sandy bottom type was assumed. 
Due to the upward refracting characteristics if the Arctic, 
the bottom type did not play a role in propagation. 

The BELLHOP MODEL, which is a Gaussian beam propagation 
model using coherent addition of ray paths, was used in 
order to keep the model run time reasonable. The 3.5 kHz 
source was placed at a depth of ten meters. Figures 32 and 
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33 are the resulting ray traces with transmission loss (TL) 
referenced dB re 1 micro Pascal. 
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Figure 32. Sound propagation through SVP with no 

staircases. 
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Figure 33. Sound Propagation model through SVP with 

staircase structure included. 


Logarithmic scaling makes it difficult to easily 
compare the two simulations, so the difference between the 
two cases is plotted to provide clarity. From Figure 34, it 
can be determined that within a range of 10000 meters, 
(range of operational detection), there is a significant 
impact from the staircases, with differences in TL up to 15 
dB. It can be noted that there is very little difference in 
the ray traces within the first 4000 meters. This is an 
expected result because in order for the staircases to have 
an impact, the sound has to be refracted so that it 
encounters the staircases at near-horizontal angles (a 
reasonable assumption for any source placement above the 
staircases because of upward refraction) The difference in 
TL using this idealized model is between -15 and 15 dB. 
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Figure 34. Difference in TL between profiles with and 

without staircases. 

Considering the idealized situation, it is reasonable 
to assume that scattering from the ice, or other ocean 
sounds speed structure, could reduce or enhance the range of 
impact of the staircases. It should be noted that even a 3dB 
difference can cause a significant impact on detection and 
should be considered. Further research is needed to realize 
that full impact of these staircases on sound propagation. 
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VI. CONCLUSIONS 


Among various aspects of the diffusive convection 
explored in this study, the most significant was the 
development of the layer-merging theory for thermohaline 
staircases. The theory based analysis of the Beaufort Gyre 
made it possible to evaluate existing flux laws for the 
ocean. From ice-tethered profiler data, there is evidence 
that merging follows Radko (2007) B-merger theory. Applying 
the theory directly to in situ data, it became evident that 
the pattern of the normalized growth-rate was the same. 
However, it was determined that the laboratory derived flux 
laws cannot be directly applied to Arctic staircases. A 
coefficient is necessary to convert from laboratory to the 
ocean. The coefficient is determined by analyzing the 
merging pattern of diffusive layers, and relating it to 
theory based flux laws. The proposed theoretical model of 
layer merging events is consistent with the field data. 

Based on this consistency, a coefficient "A" 0(1) can 
be determined and applied to better correlate laboratory 
flux laws with ocean data. Using different expressions for 
the Rp dependent coefficient of the 4/3 flux law (Kelley, 

1990; Marmorino, 1976), we estimate that the average flux 
from the thermohaline staircases is from 1-6 Wm''-2. This 
range of fluxes is significantly higher than if the original 
4/3 flux law is applied. Inconsistencies in the results of 
laboratory-derived fluxes provide significant room for error 
in direct extrapolation to the ocean. Therefore, the higher 
fluxes derived through merging theory cannot be discounted. 
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Fluxes in the Arctic are currently of great concern. As 
Turner and Veronis (2004) suggest through their laboratory 
experimentation on ice melt, double-diffusive processes may 
be somewhat responsible for melting Arctic ice. While it 
has been suggested that fluxes in the Beaufort Gyre 
staircases cannot penetrate the halocline, it is entirely 
possible that regions of double-diffusive convection, where 
the halocline is weaker, could have a greater effect. It is 
also noted that the heat budget for the Arctic has 
approximately 1/3 of its fluxes not accounted for by solar 
input (Perovich et al. 2001) . While evidence is 
inconclusive as to the exact contribution of convective 
diffusive fluxes to mixed layer over-all flux, it would be 
negligent not to consider the possibility of its 
contribution. 

In the cycle of global climate change, attention is 
once again being diverted to the Arctic. The accessibility 
of resources and shipping routes that a melted ice cap would 
provide necessitate an operational focus on the impact of 
microstructures on sound propagation. An idealized 
propagation model was used to determine whether or not 
further research was necessary, and found that under normal 
operational conditions, there was an impact from the 
staircase structure of up to 30dB. While further research 
could increase or decrease this magnitude, the difference in 
transmission loss should be considered in future Arctic 
operations, and should be further studied to determine the 
actual magnitude of impact from this structure. 


54 



LIST OF REFERENCES 


Chin-Bing, S. A., King, D. B., & Boyd, J. D. (1994) The 
Effects of Ocean Environmental Variability on 
Underwater Acoustic Propagation Forecasting. In 
Robinson, A. & Lee, D., Oceanography and Acoustics 
Prediction and Propagation Models (24-49). New 
York: American Institute of Physics. 

Crapper, P. F. (1975). Measurements across a diffusive 
interface. Deep-Sea Research, 22 , 537-545. 

Del Grosso, V. A. (1974) . New equation for the speed of 
sound in natural waters (with comparisons to other 
equations). Journal of the Acoustical Society of 
America 56 ( 4 ), 1084-1091. 

Fernando, H. J. S. (1987) . The formation of layered 

structure when a stable salinity gradient is heated 
from below. Journal of Fluid Mechanics, 182 (425-442), 
525-541. 

Huppert, H. E. (1971) . On the stability of a series of 
double-diffusive layers. Deep Sea Research, 18 , 1005- 

1021 . 

Johnson, G. C., Toole, J. M., & Larson, N. G. (2007) . 

Sensor corrections for Sea-Bird SBE-41Cp and SBE-41 
CTDs. Journal of Atmospheric and Oceanic Technology, 
24 , 1117-1130. 

Kelley, D. E. (1988). Explaining effective diffusivities 
within diffusive oceanic staircases. In J. C. J. 
Nihoul, & B.M. Jamart (Eds.), Small-scale turbulence 
and mixing in the ocean (481-502). Amsterdam: 
Elsevier. 

Kelley, D. E. (1990). Fluxes through diffusive staircases: 
a new formulation. Journal of Geophysical Research, 
95 , 3365-3371. 

Kelley, D. E., Fernando, H. J. S., Gargett, A. E., Tanny, J. 
& Ozsoy, E. (2003). The diffusive regime of double- 
diffusive convection. Progress in Oceanography. 56 , 
461-481. 


55 


Krishfield, R., & Coauthors (2006): Design and operation of 
automated Ice-Tethered Profilers for real-time seawater 
observations in the polar oceans. Woods Hole 
Oceanographic Institution Tech. Rep. WHOI-2006-NN, 30 

pp. http://hdl.handle .net/1912/1170 , accessed September 
12, 2007. 

Marmorino, G. 0. & D. R. Caldwell (1976) . Heat and salt 

transport through a diffusive thermohaline interface. 
Deep-Sea Research, 23, 59-67. 

Padman L. & Dillon, T. M. (1987) . Vertical heat fluxes 
through the Beaufort Sea thermohaline staircase. 
Journal of Geophysical Research, 92 (CIO), 10799-10806. 

Perovich, D. K., Grenfell, T. C., Richter-Menge, J. A., 
Light, B., Tucker, W. B., & Eicken, H. (2003) . Thin 

and thinner: Sea ice mass balance measurements during 

SHEBA. Journal of Geophysical Research, 108(C3), 1-21. 

Radko, T. (2005) . What determines the thickness of layers 
in a thermohaline staircase? Journal of Fluid 

Mechanics, 523, 79-98. 

Radko, T. (2007) . Mechanics of merging events for a series 
of layers in a stratified turbulent fluid. Journal of 
Fluid Mechanics, 000, 1-23. 

Stommel, H., Arons, A. B., & Blanchard, D. (1956). An 
oceanographical curiosity: the perpetual salt 
fountain. Deep-Sea Research, 3, 152-153. 

Stern, M. E. (1960). The 'salt fountain' and thermohaline 
convection. Tellus, 12, 172-175. 

Turner, J. S. (1965). The coupled turbulent transports of 
salt and heat across a sharp density interface. 
International Journal of Heat and Mass Transport, 8, 
759-767. 

Turner, J. S. & Veronis, G. (2004) . The influence of 

double-diffusive processes on the melting of ice in 
the Arctic Ocean: laboratory analogue experiments and 
their interpretation. Journal of Marine Systems, 45, 
21-37. 


56 



Turner, J. S. (1968). The behavior of a stable salinity 
gradient heated from below. Journal of Fluid 
Mechanics, 33 , 183-200. 

Turner, J. S. (1973) Buoyancy Effects in Fluids. New York 
Cambridge University Press. 


57 



THIS PAGE INTENTIONALLY LEET BLANK 


58 



INITIAL DISTRIBUTION LIST 


1. Defense Technical Information Center 
Ft. Belvoir, Virginia 

2. Dudley Knox Library 
Naval Postgraduate School 
Monterey, California 

3. Dr. Mary Batteen 

Naval Postgraduate School 
Monterey, California 

4. Dr. Timour Radko 

Naval Postgraduate School 
Monterey, California 

5. Commander Ben Reeder, United States Navy 
Naval Postgraduate School 

Monterey, California 


59 



